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^ ■ Abstract 

• In this paper we study the gauge mediation signal with the ATLAS detector at the CERN 

. LHC. We focus on the case where the NLSP is the long-hved hghtest neutrahno (x?) which decays 

dominantly into a photon (7) and a gravitino (G). A non-pointing photon from the neutrahno 

decay can be detected with good position and time resolutions by the electormagnetic calorimeter 

■ (ECAL), while the photon momentum would be precisely measured if the photon is converted inside 

O '. 

. the inner tracking detector before reaching the ECAL. A new technique is developed to determine 

o 



the masses of the slepton {€) and the neutrahno from events with a lepton and a converted non- 
pointing photon arising from the cascade decay I ixi — > i'yG- A Monte Carlo simulation at a 

Oh! 

D . sample point shows that the masses would be measured with an error of 3% for O(IOO) selected £7 



p. 

pairs. Once the sparticle masses are determined by this method, the decay time and momentum 
^ ' of the neutrahno are solved using the ECAL data and the lepton momentum only, for all £j pairs 

without the photon conversion. We estimate the sensitivity to the neutrahno lifetime for cr = 10 cm 
to 0(10) m. 

PACS numbers: 12.60.Jv, 14.80.Ly 
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I. INTRODUCTION 



Minimal Supersymmetric Standard Model (MSSM) is a promising candidate beyond the 
Standard Model. As the supersymmetry (SUSY) must be spontaneously broken, the MSSM 
needs an additional sector (the hidden sector) which breaks the supersymmetry while avoid- 
ing the FCNC problem. The origin of the SUSY breaking and the mediation to the MSSM 
sector arc therefore the key feature of SUSY models. The hidden sector SUSY breaking are 
expressed in terms of the order parameter of the SUSY breaking F and the scale of the SUSY 
breaking mediation to the MSSM sector M. The mass scale of MSSM sparticles Msusy is 
then of the order of XF/M, where A is the couphng of the hidden sector to the MSSM sector. 
If M ~ Mpi, Msusy = 1 TeV corresponds to VF ~ 10^° GeV. This class of mediation is 
called the supergravity (SUGRA) model. On the other hand, the SUSY breaking mediation 
may be due to renormalizable interactions, such as the gauge interaction. This is called 
"gauge mediation" (GM) models [1, 2]. In the GM models M and F are arbitrary and we 
expect M < Mpi. 

The GM models are described by a few parameters. The MSSM gaugino masses Mj 
(i = 1, 2, 3) and slepton masses are of the order of ctjF/M in the simplest GM model, where 
ai denotes each gauge coupling constant. On the other hand, the gravitino mass ttIq is 
proportional to Fq/M^i where Fq is the order of the SUSY breaking of the total system 
(Fo>F) [3, 4]. 

Because M < Mpi, the lightest SUSY particle (LSP) is the gravitino (G) in the GM 
models. The next hghtest SUSY particle (NLSP) is a particle in the MSSM sector which 
decays into a gravitino [5] . If the lightest neutralino (xi) is the NLSP, the dominant decay 
mode is x? iG. The neutralino lifetime ct^ is a function of Fq and m^o, and may be 
long-lived. 

The CERN large hadron collider (LHC) is a pp collider at center of mass energy of 14 TeV. 
The LHC is expected to start its physics runs in 2007. The initial integrated luminosity 
will be 10 fb~^/yr at the beginning (low luminosity runs), and then upgrade to 100 fb^^/yr 
(high luminosity runs). Signatures of the GM models at the LHC are spectacular [6-8]. In 
the case of the neutralino NLSP, SUSY events have nearly always hard photons, which may 

^ In this paper the neutrahno hfetime is multiphed by the Hght velocity c to have a dimension of length. 
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not be pointing to the interaction point (non-pointing photons). 

In this paper we propose a new approach to study the signature of the GM models using 
the ATLAS detector at the LHC for the case where the neutrahno NLSP dominantly decays 
into a photon and a gravitino with cr longer than 0(10) cm. We use two newly developed 
techniques. One is to determine the direction of the gravitino momentum by using the arrival 
position, arrival time, and momentum of the non-pointing photon, which are measured at 
the electromagnetic calorimeter (ECAL). The precision of the photon momentum would 
be significantly improved if the non-pointing photon is converted into an e+e~ pair in the 
inner tracking detector located inside the ECAL. The other technique is the "mass relation 
method" , a mass reconstruction technique which does not rely on the conventional cndpoint 
measurement [9]. In this method wc use the fact that each event from a same cascade decay 
satisfies mass shell conditions of intermediate particles. These techniques are described in 
Sec. II and a simulation is carried out in Sec. III. A fast Monte Carlo simulation shows 
that, using the new techniques, statistical error of the masses of i and x5 can be a few % 
for C(IOO) selected ^7 pairs from the decay chain £ — > ixi — > i^G. Although more works, 
especially full detector simulations, are needed to establish the techniques, the result of the 
fast simulation is quite encouraging. 

In Sec. IV we show that events with £7 pairs are fully reconstructed by using the measured 
mass and the ECAL information. This measurement is utilized to determine the neutralino 
lifetime for 10 cm < cr < 0(10) m. This analysis does not require the photon to be converted 
in the inner detector. Therefore available number of events is significantly larger. Finally 
in Sec. V we discuss how these measurements would be translated into the fundamental 
parameters Fq, F and M in the GM models. 

II. KINEMATICS OF THE EVENTS WITH NON-PROMPT DECAY 

We first discuss the decay kinematics of the neutralino with an unnegligible lifetime. In 
FIG. 1 we schematically show a neutralino decaying into a gravitino and a photon in the 
ATLAS detector. The neutrahno is produced at the interaction point O at the time t = 0, 

and then flics to the decay point D, where the neutralino decays at t = to into a photon and 
a gravitino. The photon goes to the point A in the ECAL at t = t^. Wc dcflnc two angles 
a and ip, the former between the photon momentum and the position vector — OA, 
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FIG. 1: Decay kinematics of the NLSP (the hghtest neutrahno Xi) the ATLAS detector. 

and the latter between and the gravitino momentum p^. 

We can experimentally measure a, and the distance L = \0A\. The angle ip can then 
be calculated from the three observable as 

dry — L COS a . . 

where 4 = 7 • (1) 

L sm a 

Because the three momenta p^^, and Pq are on a same plane, the direction of the gravitino 
momentum is completely determined^. 

The information of the gravitino direction may be used to determine the sparticle masses. 
We describe this idea for the decay chain £ — ^ (ixi ~^ (-"^G- The slepton C may be copiously 



^ The formula can be shown quite easily by going into the frame where the interaction point and the 
detection point is same. We define the rapidity of the gravitino and photon, taking the boost direction 
as the z direction. To obtain the formula Eq. (1) one then goes back to the laboratory frame noting the 
additivity of the rapidity. We thank Dr. Odagiri for pointing out this. 
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produced at the LHC from X2 or Xi decays, where X2 xf ^-^e dominantly produced 
from gluino or squark decays. The neutrahno and slepton masses {m^o and m^) are related 
by the following formula; 

= 2E^E^{1 - cos^) + 2EeE^{l - cos^^g) 
+ 2EeE^{l - cos 6*^^) 

= + 4(1-00.;) ) + ^^'^-(1 - 

= amlo + 6 , (2) 
where we use the relation 

= {Pj+PgT 

= 2E^Eq{1 - cos^)) , (3) 

and neglect the lepton and gravitino masses. A pair of parameters (a, b) can be calculated 
event by event from the momenta of the lepton and the photon, and the direction of the 
gravitino momentum. Because the sparticle masses m^o and should be common for all 
events, one can determine them if we have at least two tagged events. 

In Refs. [9, 10] the masses m^o and mi are determined by measuring endpoints of mass 
distributions for events containing which come from the decay chain X2 -^-^ -^-^Xi ~^ 
££^G. Endpoints in the distributions of invariant masses rriie, rrii^, and mu'y are combined 
to solve m^, m^o and itl.^o. Note that only the events near the endpoints contribute to the 
mass determination in the endpoint analysis. 

Our proposal is quite different from the endpoint analysis. We assume a set of events 
come from a common cascade decay, and use the mass shell condition of the sparticles 
involved in the cascade decay. Each event gives an independent constraint to the masses as 
given in Eq. (2), and contributes to the mass determination. We call this " the mass relation 
method". This technique may be applied for other cascade decays of SUSY particles, which 
will be discussed in future publications. 

The ATLAS detector at the LHC has a good capability to measure non-pointing photons, 
where the barrel part of the electromagnetic calorimeter (ECAL) and the transition radiation 
tracker (TRT) will play important roles (FIG. 1). The barrel ECAL is a liquid- Argon 
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calorimeter, and covers the psuedo-rapidity range I77I < 1.4. The inner radius of the ECAL 
is 150 cm. We assume the angular resolutions of the photon arrival point at the ECAL inner 
surface to be ~ 0.004 and ~ 0.002. The ECAL energy resolution is expected to be 
(Te^/E^ = where the photon energy is given in GeV. The longitudinal and 

transverse segmentation of the ECAL gives a measure of the development of electromagnetic 
showers. The first longitudinal sampling is finely segmented in the r} direction, resulting in 
a good angular resolution of — 60 mrad/ ^/E^^, where 9 is the polar angle of the photon 
momentum with respect to the beam axis. The azimuthal angle (p of the photon momentum 
is only poorly measured by the ECAL, as the segmentation is very coarse in this direction. 
The ECAL also has an excellent time resolution, < 100 ps for > 30 GeV, confirmed 
by a test-beam experiment. 

When a photon is pointing to the interaction point, the photon momentum is precisely 
determined by the ECAL only, namely by measuring the energy deposit and the arrival 
position. However, in the case of the CM models, the photon is in general non-pointing, and 
the transverse components of the photon momentum are only poorly measured. Fortunately, 
the barrel TRT is located inside the barrel ECAL as a component of the inner tracking 
detector. This detector covers the radial range from 56 to 107 cm and the pseudo-rapidity 
range \r}\ < 0.7. As the straw tube trackers of the barrel TRT are parallel to the beam 
axis, trajectories of charged particles are precisely measured in the r-0 plane. If a photon 
is converted into an e+e~ pair before escaping the barrel TRT, the angle of the photon 
momentum can be very precisely measured. As the angle resolution is much better than the 



6 angle resolution by the ECAL, the resolution of the angle a becomes cxa = y cr^ + cr^ ~ ae. 
The material thickness of the TRT is about 10% of one radiation length at ?] ~ 0. Namely 
~10% of photons will be converted in the TRT. 

III. RECONSTRUCTION OF GRAVITINO DIRECTION AND SPARTICLE 
MASSES 

In order to test the new techniques we perform a Monte Carlo event simulation at CM 
point Gl [9, 10] with Fq, or equivalently the neutralino lifetime cr, being the only free 
parameter. The model parameters and some of sparticlc masses arc listed in TABLE I. The 
low energy SUSY parameters are calculated by ISASUSY 7.51, and the mass spectrum, the 




6 



TABLE I: Model parameters and some of the sparticle masses at point Gl. The parameter N is 
an integer number which appears in Eq. (10). 



Parameters Sparticle masses (GeV) 

F/M = 90 TeV m-g = 720 m^-^ = 324 m^o = 117 

M = 500 TeV mgr = 958 m„- = 162 m^o = 217 

A2 

= 1 m-qj, = 915 mf^ = 831 m^o = 420 

tan /? = 5 m^^ = 909 m^o = 442 
At > 



couplings and the decay branching ratios are interfaced to HERWIG version 6.4. We generate 
10^ SUSY events at this point, corresponding to an integrated luminosity of 13.9 fb~^. When 
we simulate the events, we keep the lightest neutralinos (NLSPs) stable at the generator level. 
Then a fast detector simulator ATLFAST is used for all particles except the neutralinos. The 
decay of the neutralinos and the photon conversions are simulated at the analysis stage. The 
photon conversion probability is estimated based on the detector thickness of the TRT [9] . 
If (a) a neutralino decays into a photon and a gravitino before escaping the TRT region, 
(b) the photon points to the barrel ECAL, and (c) the photon is converted inside the barrel 
TRT, then the energy, position, time, and direction of the photon are smeared by the 
Gaussian distribution according to their resolutions. The detector resolutions assumed in 
the simulation are hsted in TABLE II. We assume the time resolution of the ECAL to be 
constant — 100 ps for > 30 GeV. 

For a moment we set the neutralino lifetime to be cr = 100 cm. In this case the neutralinos 
efficiently decay in the TRT, as the outer radius of the TRT is roughly 100 cm. We first 
apply pre-selection cuts to suppress the Standard Model background: 

i) Meff > 400 GeV, 

ii) E^''' > O.lMeff. (4) 

The missing transverse energy E^^^^ is calculated from the reconstructed jets, leptons, pho- 
tons and unreconstructed calorimeter energies. The effective mass is defined by the sum of 
the missing transverse energy and the transverse momenta of the four hardest jets: 

Meff = + Pt,1 + PT,2 + PT,3 + PT,4 (5) 
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TABLE II: Detector resolutions assumed in our Monte Carlo simulation of non-pointing photons. 
The photon energy is given in GeV. As for the (p angle of the photon momentum, we assume 
the photon conversion in the TRT detector. 



Observable 




Detector 


Resolution 


Photon energy 


E-y 


ECAL 




Photon arrival time 




ECAL 


100 ps 


Photon arrival position 


V 


ECAL 


0.002 




o 


ECAL 


0.001 


Photon momentum 


e 


ECAL 


0.060/y/E^ 




4> 


TRT 


0.001 



We do not include photons from the neutralino decays in the above calculation. The effi- 
ciency of the pre-selection cuts for the generated SUSY events is 80%. 

The following cuts^ are then applied to select good 'non-pointing' photons with conversion 
in the TRT: 

1) > 30 GeV, 

2) a > 0.2, 

3) Ai^(= - L/c) > 1.0 ns. (6) 

The distributions of E'-y, a and At^ are shown in FIGs. 2(a)-(c), where the cuts l)-3) are 
sequentially applied. 

In FIG. 2(d) we plot A'^ = ip — ■^truej where ip is calculated from the measured L, a, and 
t-y using Eq. (1) and ■0true is the true value obtained from the generator information. The 
resolution is better than 40 mrad in this case. 

In order to determine the masses of the slepton and the neutralino which appear in 
the cascade decay £ — > ixi — ^ ^7^, isolated leptons (electrons and muons) with transverse 
momentum larger than 20 GeV are searched for to make a pair with the non-pointing photon. 
If there are several leptons in an event, we choose the ^7 pair which minimizes the invariant 
mass m^-y. The parameters a and b are calculated for each ^7 pair. The scatter plot in the 

^ Isolation from tracks/clusters is yet to be examined. 
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FIG. 2: Distributions of (a) Ej, (b) a, and (c) At^ at point Gl with cr = 100 cm, where the cuts 
l)-3) in Eq. (6) are sequentially applied, (d) Distribution of Aip = ijj — V'true after applying the 
cuts, where i/j is the angle between the momenta of the photon and the gravitino. The result of a 
Gaussian fit is also shown. 

(a, b) plane is shown in FIG. 3(a), where the sample contains about 120 £7 pairs. The points 
are clearly concentrated along a line. The section of the b axis and the slope of the line must 
correspond to m| and m|o because of the relation b — m| — mioa . In FIG. 3(b) we take a 
simple average of (a, b) for the events between the two dotted lines in FIG. 3(a) by dividing 
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FIG. 3: (a) Distribution of ^7 pairs in (a, b) plane. The region between dotted lines shows our 
selection cut. (b) Points show the average values of (a, b) between the two dotted lines in (a). The 
region is divided into 9 bins. The solid line shows a result of a linear fit. 

the region of a into 9 bins, and we fit the averaged data by the hnear function Eq. (2). The 
fit results are = 162.1 GeV and m^o — 117.5 GeV, while the input values are 161.7 GeV 
and 117.0 GeV, respectively. 

To estimate errors of the mass measurement, the simulation and the fit are repeated 100 
times with different random number seeds. The fitted masses are plotted in FIG. 4. By fitting 
the distribution, wc obtain the errors of the masses as a^- = 2.7 GcV and Cm,,, = 3.5 GeV. 
They correspond to relative mass errors of 2 ~ 3%. If one of the sparticle masses is precisely 
determined by some other measurements, the other sparticle mass can be determined with 
a precision of ~ 300 MeV by the mass relation method. 

IV. FULL RECONSTRUCTION AND LIFETIME MEASUREMENT 

In this section wc demonstrate full reconstruction of the cascade decay i £xi ~^ ^iG, 
and we show that the neutralino lifetime can be determined by the reconstructed decay 
time and momentum of the neutralino. This analysis becomes possible after the precise 
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FIG. 4: Distributions of the fit results of (a) the slepton mass and (b) the neutralino mass rn^o. 
The simulation and the fit are repeated 100 times with different random number seeds. Results of 
Gaussian fitting are also shown. 

determination of and m^o described in the previous section. 

Here we study events with leptons and non-pointing photons, where the photons may or 
may not be converted in the inner detector. For each i'j pair, the arrival time t^, arrival 
point x~f and energy E.y of the photon, the longitudinal component of the photon momentum, 
and the lepton momentum pi would be directly measured, while the transverse components 
of the photon momentum, the neutralino four momentum and its decay time t^, (< t^) 
are not directly measured. On the other hand, we have the following equations involving 
the unknown quantities. 



where and Uy are velocity vectors (|uy| = c), and the lepton and gravitino masses are 
neglected. Provided that we know the two masses mi and m^o, we can solve all the unknown 





(7) 
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parameters from the equations. There are two solutions for each £7 pair, which we obtain 
by numerically solving Eq. (7). 

The simulation and analysis are modified from those of the previous section. The neu- 
tralino must decay inside the ECAL and the photon from the neutralino decay must enter 
the barrel ECAL, resulting in a much larger fiducial decay volume. We assume the fiducial 
volume as r < 150 cm and \z\ < 300 cm, where r and z are the radial distance from the 
beam axis and the distance from the interaction point along the beam axis, respectively. In 
addition, as the photon conversion is not required in the inner detector, the acceptance of 
the events further increases. After applying the pre-selection cuts given in Eq. (4), we select 
non-pointing photons by the following cuts; 

1) E^>30 GeV, 

2) > 0.2 rad, 

3) At^ > 1 ns. (8) 

Here = 9^^ — is the difference between the polar angles of and x^, where the 
polar angles are measured from the beam axis. Note that the cut 2) is different from that 
in Eq. (6), because one cannot determine the angle a without precise measurement of the 
(j) angle using the photon conversion. This A^ cut is less efficient than the a cut in Eq. (6). 
However, the overall selection efficiency of this analysis is much larger because we have the 
larger fiducial decay volume and do not require the photon conversion in the inner detector. 

Isolated leptons with transverse momentum larger than 20 GeV are used to make £7 
pairs. For each i'j pair, there are two solutions for decay kinematics. We take solutions if 
the reconstructed decay points are in the fiducial volume (e.g. r < 100 cm and \z\ < 300 cm). 
Therefore each £7 pair may be counted twice. If there are several leptons in an event, all 
leptons are tried for each non-pointing photon, and accepted if the ^7 pair satisfies the 
fiducial condition. 

In FIG. 5(a) we show a scatter plot of the reconstructed decay time to and the true 
decay time for cr = 100 cm. Here we use the input and m^o for reconstruction. This 
plot shows that the neutralino decay time is correctly reconstructed for a significant fraction 
of the events. We also show the distribution of ^£1/7^^ in FIG. 5(b) for cr = 50 cm, 100 cm, 
and 200 cm, where the factor 1/7^ corrects the effect of Lorenz boost. The distribution shows 
the expected exponential damping toward large ^d/Tx values. The geometric acceptance of 
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FIG. 5: (a) The distribution of the neutrahno decay time t£) vs the true value, (b) Distributions 
oi to/jx ~ ^0 (dashed), 100 cm (sohd) and 200 cm (dotted). 

the events would be small when cr7^ > i?j„ or t'J^ < 1 ns, where (~ 150 cm) is the 
inner radius of the ECAL and 1 ns is our cut value on the difference of the photon arrival 
time. 

The geometric effect to the acceptance may be corrected by a study with full detector sim- 
ulations. The momentum distribution of the neutralino would be an important uncertainty, 
as the transverse momentum distribution should depend on the gluino and squark masses. 
The heavier sparticle masses may be measured precisely by the jxi or jjXi invariant mass 
distribution for selected jets (j) and a neutralino, where the neutralino momentum is recon- 
structed by solving Eqs. (7). Therefore, we assume the systematics due to the momentum 
distribution of the neutralino would be small enough. 

Assuming that the systematic errors are controlled, we may study the sensitivity to the 
lifetime cr using the two measured values; 

(a) (^d/Tx)- the average of the corrected decay time, 

(b) NiY- the number of accepted £j pairs. Note that we do not count an ij pair more 
than once to avoid over-counting. 
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FIG. 6: (a) Average {toMxl ^) ^ functions of cr for an integrated luminosity of 13.9 
fb~^ at point Gl. A dot and an error bar show the mean value and the standard deviation of 100 
simulations at each cr value, respectively. 

They are plotted as functions of cr in FIG. 6, where we repeat the simulation for an inte- 
grated luminosity of 13.9 fb~^ hundred times to obtain the mean (to/jx) values 
and their standard deviations. In the plot, (io/Tx) larger than 1 ns even for cr = 10 cm 
because of the cut Ai^ > 1 ns. Note that we do not have sensitivity on ^0 when cr < 10 cm 
because the time resolution is ~0.1 ns, and c x 0.1 ns ~ 3 cm. The number of £7 pairs with 
photon conversion (A^™"^^) is 13.0 and 80.5 for cr = 10 cm and cr = 30 cm, respectively. As 
the number of converted photons is rather small for cr = 10 cm, a large integrated lumi- 
nosity is needed to determine the sparticle masses. For large cr values, the average (to/jx) 
is saturated since most of the neutralinos decay outside the detector. Indeed the number of 
the reconstructed events takes its maximum at cr ~ 100 cm and decreases monotonically as 
shown in FIG. 6(b). 

The sensitivity of the {tD/lx) measurement to the lifetime cr is estimated by the error of 
the measurement ^{to/^x)- Namely we define the error of the lifetime Acr by the following 
formula: 

g{cT ± Act) = {tohx) ± HtD/lx) , (9) 
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FIG. 7: Estimated resolution of the lifetime cr for an integrated luminosity of 13.9 fb^^ from (a) 
the average Id/'^x (^) number of ^7 pairs N^^. The input and rriy. are used for the 
reconstruction and their errors are ignored. 

where g{cT) is a function to describe the average {to/lx)- "^^^^ function is numerically 
obtained by fitting the {to/jx) ^ second power polynomial function of cr using the 
average values in the region within cr ± 50 cm. When cr ^ 100 cm, most of the neutralinos 
decay after reaching the ECAL and (^d/Tx) looses the sensitivity to cr. On the other hand, 
the number of the £7 pairs N^^ sensitively reduces when cr is longer than 100 cm. We can 
estimate the sensitivity of N^^^ to cr in a same way. The result is shown in FIG .7 and 
summarized in TABLE III. Note that we ignore the errors of and m^o for a moment 
and use the input mass values. Note also that only statistical errors are considered here. 
Systematic errors are yet to be studied with full detector simulations. 

For large cr values the number of selected events is limited. In TABLE IV we summarize 
Ne^ and N^^""" for an integrated luminosity of 13.9 fb^^ For cr < 2000 cm the number 
iV|°'^^ still exceeds 130 events if we assume an integrated luminosity of 100 fb^^ at point Gl. 
Therefore, it might be still possible to determine the neutralino and slepton masses. 
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TABLE III: Relative resolution of cr. As in FIG. 7 (a) and (b) are estimated from the average 
{tD/lx) number of ^7 pairs N^(^, respectively. 



(a) Prom {to/lx) (b) I^om N^^ 



(cm) 


/^ct/ct 


cr (cm) 




10 


0.173 


200 


0.062 


20 


0.062 


250 


0.038 


30 


0.041 


300 


0.034 


40 


0.036 


350 


0.036 


50 


0.037 


400 


0.030 


60 


0.027 


450 


0.031 


70 


0.037 


500 


0.040 


80 


0.039 


600 


0.025 


90 


0.043 


700 


0.037 


100 


0.045 


800 


0.037 


150 


0.043 


900 


0.037 


200 


0.047 


1000 


0.039 


250 


0.061 


1500 


0.048 


300 


0.070 


2000 


0.033 



TABLE IV: Number of ^7 pairs for an integrated luminosity of 13.9 fb ^ . is the number of 

]")Mirs with photon is couA-orsion in IIk^ irmer detector. 



cr (cm) A^^°"^ 

10 147 13.0 

30 1424 80.5 

100 3193 116.3 

300 2413 72.8 

2000 536 14.3 
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V. DETERMINATION OF THE FUNDAMENTAL PARAMETERS 



In this section we utilize the measurement of the sparticle masses and the neutrahno 
hfetime to determine more fundamental parameters in the GM scenario. Prom the measure- 
ment of the slepton and neutralino masses, the ratio F/M is determined with a precision of 
a few %. This can be seen in the expression of the masses at the messenger scale as [11]; 



mi = 



3aj{M) 

407r2 M2 

m? 



N X f [F/M^) 



? = + MM^ X / (FIM^) (10) 

where N is an integer number, while g and / are some functions which satisfy /(O) = g{Q) ~ 1 
for F <^ M^. The sparticle masses are proportional to F/M in this limit. 

The absolute size of M (or F) is rather difficult to determine because it only appears 
through the sfermion mass running from M to the SUSY scale. A study [9] shows that the 
relative error AM/M is ~ 30% with 1% mass errors at this point. 

The neutralino lifetime depends on the order parameter of the total SUSY breaking Fq 
and the neutralino mass m^o [5]; 

100GeV\V a/F 



^7 V ^xS ' ^100 TeV ' ' ^ ' 

where k.y — jA^nCOS^w^ + Ni2sm.9w\'^ with 9w being the Weinberg angle, and Nij is the 
neutrahno mixing angles. The constant is = cos^ dw for the bino-like The parameter 
Fo is also related to the gravitino mass itself [3, 4]; 



mn = ^T7^ = Z r^.r 2.4 eV. (12) 

The sensitivity to the neutrahno hfetime is about 4% for cr ~ 100 cm as given in the 
previous section. Ignoring the errors of the neutralino mass and A;^,, this is translated to an 
error of about 1% for ^/Fq. The lifetime cr = 100 cm corresponds to \/Tq ~ 1000 TeV and 
rriQ ~ 0.2 keV, respectively. 

If CT is 0(10) ~ 0(100) m, the ratio of the number of £7 pairs to that of SUGRA-hke 
events is a more sensitive measure for the lifetime. For cr = 2000 cm {-s/Fq ~ 2400 TeV 
and rriQ ~ 1.4 keV), the estimated error is about 4%. This corresponds to a 1(2)% error 
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on ^/Fq {itiq). The mass error of x5 is estimated to be around 3%, because the number 
of events with a converted photon is around 100 for an integrated luminosity of 100 fb~^, 
which is obtained in one year at the high luminosity runs of the LHC. Prom Eqs. (10) and 
(12) the error of Fq {mo) due to the mass error is 3(6)% in this case. For cr = 10^ cm 
the number of events with a lepton and a converted photon is about 70 for the ultimate 
integrated luminosity of 300 fb~^. 

The above estimate is rather optimistic because we do not consider systematic errors of 
the measurement. In addition, we have so far not considered the effect of background events 
from the Standard Model processes and the SUSY production itself, where prompt photons 
from the interaction point mimic non-pointing photons due to the limited detector resolution. 
Note that in our Monte Carlo sample, 28757 slcptons are produced, and ~ 2500 (~70) of 
them are accepted as the sample with a (converted) non-pointing photon for cr = 300 cm. 
Especially for a large cr case, the number of real non-pointing photons is small and the 
analysis might severely suffer from the background. Studies with full detector simulations are 
necessary to understand the systematic errors and the background, and should be completed 
before the start of the LHC physics run. If events with a non-pointing photon converting 
in the TRT cannot be used for mass reconstruction due to the backgrounds, the Am^o ~ 
30 % is expected from the endpoint analysis involving jets and leptons. In this case the 
uncertainty of the neutralino mass would lead to an error of the gravitino mass by a factor 
of two. 
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